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1.  Introduction 


When  first  developed  in  World  War  II,  radar  was  used  to  detect  the  presence 
of  objects  within  its  range.  Initially,  radars  could  report  the  position  and  possibly 
the  velocity  of  detected  objects  (called  “targets”),  but  it  soon  became  evident  that 
more  information  was  desirable.  Modern  civilian  air  control  solves  the  information 
problem  by  requiring  each  aircraft  to  transmit  information  about  itself,  using  a 
transponder  [1].  However,  the  problem  of  obtaining  a  characterization  of  the  target 
when  it  is  unable  or  unwilling  to  volunteer  that  information  remains  a  topic  of 
current  research. 

The  response  of  the  target  to  the  radar  signal  (the  radar  return)  contains  a 
considerable  amount  of  information  about  the  target,  because  the  target  distorts 
the  radar  waves  as  they  wash  over  it.  The  exact  nature  of  the  distortion  depends  on 
the  shape  and  orientation  of  the  target.  It  should  therefore  be  possible  to  extract 
information  about  the  shape  and  orientation  of  the  target  from  the  radar  return. 
If  we  have  a  database  containing  shape  information  about  all  known  aircraft  (the 
catalog),  we  can  identify  the  target  by  comparing  the  shape  information  extracted 
from  the  radar  return  to  the  various  entries  in  the  catalog.  This  process  is  known 
as  radar  target  identification  (RTI). 

The  radar  target  identification  system  typically  consists  of  a  signal  processing 
stage  followed  by  a  feature  classification  stage.  The  signal  processing  step  involves 
operating  on  the  raw  radar  data  to  extract  salient  “features”  of  that  data  which  can 
be  readily  used  for  target  classification.  This  report  addresses  the  signal  processing 
part  of  the  RTI  problem. 


Traditionally,  target  identification  methods  have  used  direct  frequency  domain 
data  from  the  radar.  The  “features”  used  for  classification  consist  of  raw  frequency 
domain  data  [2, 3, 4, 5].  This  method  requires  little  or  no  signal  processing,  but  has 
the  disadvantage  that  one  obtains  no  geometric  characterization  of  the  target.  (Ge¬ 
ometric  information  is  useful  because  it  can  provide  information  about  the  target 
even  when  classification  or  identification  is  not  possible.)  More  recently,  Fourier 
Transform  (FT)  methods  have  been  applied  [6];  these  methods  first  convert  the  fre¬ 
quency  domain  radar  data  into  an  estimate  of  the  downrange  impulse  response  of 
the  target.  Classification  then  proceeds  based  on  some  extracted  features  from  this 
downrange  impulse  response  (such  as  range  of  strong  scattering  centers,  etc.).  While 
this  concept  is  potentially  very  useful  for  RTI  applications,  the  use  of  FT  methods 
presents  some  problems.  First,  FT  methods  are  resolution  limited,  so  closely  spaced 
scattering  centers  may  not  be  resolved.  Second,  FT  methods  are  nonparametric; 
thus,  there  is  no  data  reduction  from  the  raw  radar  data  to  the  downrange  im¬ 
pulse  response.  Data  reduction  is  desirable  because  classification  algorithms  are 
computationally  burdensome  if  the  number  of  features  is  large. 

In  this  report  we  develop  and  test  an  alternate  radar  signal  processing  approach. 
This  approach  applies  AutoRegressive  Moving  Average  (ARMA)  modeling  tech¬ 
niques  to  radar  target  data.  Like  FT  methods,  the  ARMA  modeling  technique 
produces  an  estimate  of  the  downrange  impulse  response  of  the  target.  However, 
unlike  the  FT  model  the  ARMA  method  is  parametric;  the  output  of  the  signal  pro¬ 
cessing  stage  consists  of  a  small  number  of  parameters  which  can  be  directly  used 
for  classification.  ARMA  methods  are  also  not  resolution  limited  by  the  bandwidth 
of  the  radar  data  as  axe  FT  methods.  Thus,  these  ARMA  modeling  methods  posses 


some  potentially  useful  properties  for  radar  target  identification. 

In  this  report  we  first  derive  the  ARM  A  signal  processing  method.  We  then 
present  an  analytical  comparison  between  the  ARMA  method  and  the  FT  method. 
We  also  discuss  some  strategies  of  data  averaging  which  can  be  used  in  both  meth¬ 
ods.  Finally,  we  apply  these  signal  processing  techniques  to  compact  range  mea¬ 
surements  of  scale  models  of  several  commercial  aircraft. 
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2.  The  ARMA  Modeling  Method 


This  chapter  describes  the  basic  theory  behind  the  ARMA  modeling  method  as 
applied  to  the  RTI  problem.  Note  that  the  theory  as  presented  here  is  somewhat 
different  from  standard  developments  [7,8],  because  the  standard  development  as¬ 
sumes  that  the  measured  data  is  in  the  time  domain,  whereas  the  radar  data  is  in  the 
frequency  domain.  Thus,  in  the  radar  target  identification  application,  the  time  and 
frequency  domains  are  interchanged  with  respect  to  standard  ARMA  applications. 

2.1  The  Radar  Target  Data 

It  is  assumed  that  measurements  of  the  response  of  the  targets  are  in  the  form 
of  coherent  steady-state  response  values  (amplitude  and  phase)  at  a  number  of 
frequencies;  this  corresponds  to  frequency  measurements  taken  with  a  radar  that 
emits  either  a  continuous  wave  or  a  relatively  long  pulse  (with  respect  to  the  length 
of  the  target).  Realistic  data  for  this  type  of  measurement  is  available  from  the  OSU 
ESL  Compact  Range  Database.  Due  to  the  construction  of  the  Compact  Range, 
this  data  can  be  considered  noise  and  distortion  free,  and  is  therefore  suitable  for 
the  construction  of  a  catalog  [2].  To  simulate  the  measurements  a  real  radar  would 
take,  the  compact  range  data  can  be  corrupted  by  noise. 

The  data  from  the  Compact  Range  is  organized  as  follows:  for  each  target  there 
are  two  arrays,  amp  and  phase ,  such  that: 


amp{ 

=  \y(fi)\ 

(2.1) 

phasei 

=  Mfi) 

(2.2) 

where  y(fi)  is  the  steady-state  response  of  the  target  at  frequency  f,  =  fo  +  i6/.  For 
the  particular  data  base  we  have  used  in  this  study, 


fo  =  1GHz 
6j  =  50MHz 


(2.3) 

(2.4) 


while  i  ranges  from  0-220,  so  the  highest  available  frequency  is  12  GKz. 

The  targets  are  scale  models  of  existing  commercial  aircraft  (Boeing  747,  727, 
and  707;  DC10  and  Concorde);  scale  factors  vary  from  130  to  200,  so  the  compact 
range  measurement  frequencies  correspond  to  “real-world”  frequencies  in  the  5-92 
MHz  range). 

2.2  The  ARMA  Model 

ARMA  modeling  techniques  are  used  in  various  areas:  in  communications,  in¬ 
cluding  adaptive  matched  filters,  in  transmission  channels  (see  [9]),  and  in  geologic 
surveying  by  means  of  acoustical  procedures  (used  e.g.  to  find  oil  and  gas:  see 
[10,11,12]).  Throughout  this  report  we  will  refer  to  the  frequency  domain  represen¬ 
tation  of  the  radar  return  as  the  data ,  and  to  the  range  or  time  domain  represen¬ 
tation  as  the  downrangt  signal,  in  keeping  with  standard  terminology.  This  report 
will  be  solely  concerned  with  the  application  of  ARMA  techniques  to  the  RTI  prob¬ 
lem,  and  therefore  use  the  “reverse”  domain  formulation  that  is  appropriate  to  that 
problem;  so  data  is  transformed  from  the  frequency  domain  into  the  time  domain, 
not  the  other  way  around. 

The  signal  processing  problem  considered  here  is:  given  N  coherent  radar  mea¬ 
surements  {y(fk)}k=oi  estimate  the  downrange  impulse  response  of  the  target.  The 
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general  formulation  of  such  a  signal  estimator  entails  two  steps 

1.  Parameterize  the  impulse  response  function  (choose  a  model). 

2.  Estimate  the  parameters  of  that  model. 


The  parameterization  we  use  is  based  on  the  ARMA  model.  The  ARMA  model 


assumes  that  the  radar  data  can  be  described  as: 

na  nc 

y(fk)  =  -  +  Y1 

1=1  7=0 

where: 

fk  =  fo  +  k6f  (2.5) 

y(fk)  =  coherent  radar  data  at  frequency  fk  (2.6) 

£(fk)  =  a  complex  zero  mean  white  noise  process  with  variance  1 

(orthonormal  white  noise)  (2.7) 

This  model  can  be  represented  as: 

A(q-')v(fk)  «  C(<T’)f(A)  (2.8) 

where 

A(q  x)  =  1  +  a\q  1  +  ■  •  •  +  o,naq  na  (2.9) 

C{q~l )  =  Co  +  C\q~x  +  •  •  •  +  cncq  nc  (2.10) 

q~x  =  unit  delay  operator  (  q~ly{fk)  =  y(fk-i)  )•  (2.11) 


Note  that  we  need  to  specify  the  orders  na  and  nc  of  the  polynomials  A  and  C.  We 
will  discuss  ways  to  do  this  for  our  given  application  in  Chapter  4. 
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Prom  this  description  of  the  discrete  frequency  domain  measurements,  we  can 
find  a  description  of  the  downrange  response  via  the  inverse  2-transform: 

C(e») 


Y(  0  = 


A(ejt) 


— tt  <  l  <  n 


(2.12) 


where  /  is  scaled  time  and  is  related  to  actual  time  by 

—  /  +  7T 


t  = 


4n6t 


and  where  6j  is  the  separation  between  the  frequency  samples.  The  range  in  meters 
relates  to  the  time  delay  by  the  speed  of  light;  thus  the  range  relates  to  scaled  time 
by: 

—  /  +  7T 


r  —  c- 


47 r6f 


(2.13) 


From  here  we  can  also  very  easily  find  the  maximum  unambiguous  range:  this  is 
the  value  of  r  for  /  =  — tt: 


R  2  8S 


(2.14) 


If  the  target  extends  beyond  this  range,  aliasing  effects  will  occur.  That  is,  responses 
in  this  range  will  be  “folded  over”  and  added  to  the  response  in  the  unambiguous 
range  [13]. 

Combining  equations  (2.13)  and  (2.14),  and  substituting  into  equation  (2.12)  we 
get 

C(eJ>(1-2r/*)) 


Y ^  4(e;*(i-2r/K)) 

as  the  target  impulse  response,  where  r  is  the  range  in  meters,  or 
C(eMi-*/T)) 


(2.15) 


Y(t)  = 


A(eJ'(i-2*/T)) 


(2.16) 
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where  t  is  the  time  in  seconds,  and  T  —  Rfc  is  the  maximum  unambiguous  time. 


It  should  be  noted  from  (2.15)  that  the  estimated  downrange  response  is  a  ratio 
of  two  polynomials.  This  type  of  model  offers  the  combined  advantages  of  parameter 
parsimony  and  good  approximation  performance  [14].  In  fact,  it  can  be  shown  [14] 
that  any  function  can  be  approximated  arbitrarily  closely  with  the  model  (2.15). 


2.3  Pole- Residue  Parameterization  of  the  ARM  A  Model 


While  the  ARMA  parameters  relate  directly  to  the  shape  of  the  signal,  there 


are  other,  equivalent  parameterizations  that  are  more  directly  related  to  the  char¬ 


acteristics  of  a  radar  signal.  Specifically,  the  impulse  response  model  can  be  written 


F(r)  ~  _  p. 


(2.17) 


where  the  p;  are  the  roots  of  A(z),  and  the  d,  are  the  partial  fraction  expansion 
coefficients.  The  representation  (2.17)  is  equivalent  to  (2.15)  if  nc  =  na  —  1;  the 


only  difference  between  the  representations  is  in  the  parameters  used  to  represent 


the  model. 


The  advantage  of  representation  (2.17)  over  (2.15)  is  that  each  pole  pt  directly 


corresponds  to  a  scattering  center,  and  the  residue  relates  to  the  amplitude  as¬ 


sociated  with  this  scattering  center.  Thus,  this  parameterization  is  more  directly 


related  to  the  underlying  physical  properties  of  the  target  than  (2.15)  is. 


Each  “pole”  pi  corresponds  to  a  scattering  center  on  the  target.  The  argument 


of  pi  (argpi)  relates  to  the  range  r,-  of  this  scattering  center  by 


r.  =  ff(l  —  2—j~ ) 
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It  can  be  seen  from  (2.17)  that  r,  is  the  range  at  which  the  component  of 
the  impulse  response  achieves  maximum  amplitude.  The  magnitude  of  the  pole 
Pi  relates  to  the  distribution  in  range  of  the  energy  received  from  this  scattering 
center;  as  |p,|  — ►  1,  the  scattered  energy  becomes  more  tightly  concentrated  at  range 
r,-;  this  corresponds  to  an  ideal  point  scatterer.  For  |p,|  ^  1,  the  energy  is  spread 
over  some  range  centered  at  r,.  The  d,  parameter  gives  the  amplitude  of  the 
scattering  center  return. 

The  above  paragraph  relates  the  structure  of  the  ARMA  model  to  properties  of 
the  impulse  response  of  the  target,  i.e.  properties  in  the  time  domain.  It  is  useful 
to  consider  the  properties  of  the  model  in  the  frequency  domain  as  well,  since  these 
properties  represent  implicit  assumptions  on  the  radar  data.  The  inverse  Fourier 
transform  of  (2.17),  yields: 


na 


y(fk )  =  Y,diPki 

s=i 


(2.18) 


In  other  words,  the  ARMA  representation  assumes  that  the  radar  data  can  be  mod¬ 
eled  as  a  number  (na)  of  damped  exponentials.  Each  exponential  term  corresponds 
to  a  scattering  center.  The  energy  P,  associated  with  the  zth  scatterer  can  be  found 
from  (2.18) 


Pi  = 


Ed.<r.(p,p-)k  -  Wf 


k=0 


l  -  |p,f 1 " 
i  -  Ip.I2 


(2.19) 


2.4  Formulation  of  a  Pole-Residue  Estimator 


The  pole-residue  model  formulated  in  the  previous  section  can  be  computed 
directly  from  the  discrete  frequency  measurements,  as  detailed  below.  This  sec- 
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tion  gives  the  formulae  for  estimating  poles  and  residues  from  a  given  data  set 

{»(/*)>£?■ 

In  order  to  find  an  estimate  for  F(r),  we  need  to  find  the  A  and  C  polynomials 
or,  equivalently,  the  p,  and  d,  coefficients.  Below  we  present  a  method  adapted  from 
the  time  series  analysis  literature  [15,16]. 

The  first  step  involves  estimating  the  a,  coefficients  and  finding  the  roots  of 
the  resulting  A(z)  polynomial.  First  a  standard  estimate  for  the  autocorrelation 
sequence  corresponding  to  the  data  is  found;  here  we  use  the  standard  unbiased 


autocorrelation  estimates: 

fk=  jrzi  li  v(f, )'«(/'+>)  *  =  0,1 . it 

r-k  =  K 


(2.20) 

(2.21) 


Next,  the  a  parameters  are  estimated  from  the  r  parameters  by  solving  the  well- 
known  over  determined  Yule  Walker  equation  [17]: 


,  K  >  na  +  nc 


(2.22) 


L  rK- X  •  •  •  rK-na  J  [  Gna  J  |_  TK  J 

Now  we  can  find  the  roots  of  A(z)  using  standard  complex  polynomial  root 

finding  techniques. 

Once  the  poles  are  obtained,  the  amplitudes  d,  can  be  estimated  using  a  least 
squares  technique.  From  equation  (2.18),  such  an  estimator  is  given  by  minimizing 
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(2.23) 


or 


Pd  —  y  =  e 

This  gives  the  formula 
d  =  (P*P)~1P*y 


(2.24) 


where  P*  is  the  complex  conjugate  transpose  of  the  matrix  P.  Equations  (2.20), 
(2.22),  and  (2.24)  comprise  the  ARMA  estimation  procedure. 


2.5  Comparison  of  the  ARMA  and  DFT  Methods 


Both  the  Discrete  Fourier  Transform  and  the  ARMA  modeling  methods  perform 
the  task  of  transforming  frequency  domain  radar  data  into  a  downrange  impulse 
response.  In  this  section  we  compare  and  contrast  the  two  methods. 

The  finite  Fourier  Transform  is  given  by  [8] 


N- 1 


Wr)  =  E  »(A)ei'(‘‘!*) 

fc= 0 


Note  that  the  finite  FT  is  a  continuous  function  of  r  just  as  the  ARMA  repre¬ 
sentation  is.  The  Fast  Fourier  Transform  is  a  computationally  efficient  method  of 
computing  K(r)  at  the  N  equally  spaced  points: 


n 


rn  = 


"  (N  ~  1) 


R 
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from  the  N  frequency  measurements.  If  more  points  axe  desired,  the  data  can  be 
zero  padded;  in  the  limit  of  infinite  zero  padding,  the  FFT  gives  the  finite  Fourier 
Transform.  The  ARMA  method  gives  an  equation  for  Y(r ),  which  we  can  then 
evaluate  at  as  many  points  as  we  need. 

One  advantage  of  the  ARMA  method  is  that  it  is  not  limited  in  resolution 
by  the  bandwidth  of  the  frequency  measurements,  as  the  FT  is  [8];  that  is,  the 
ARMA  methods  are  capable  of  resolving  two  scatterers  no  matter  how  close  they 
axe  together,  irrespective  of  the  bandwidth  of  the  measurements.  The  resolution 
limit  of  the  DFT  is  caused  by  the  implicit  assumption  of  the  DFT  that  the  data 
is  identically  zero  outside  the  measurement  interval.  The  ARMA  method,  on  the 
other  hand,  implicitly  assumes  that  the  data  outside  the  measurement  interval  is  a 
continuation  of  the  data  in  that  interval,  that  is  the  equation  (2.18)  is  assumed  to 
hold  for  all  k,  not  just  0  <  k  <  N  —  1  (see  chapter  3  for  details). 

Another  advantage  of  the  ARMA  methods  is  that  it  produces  the  estimates  in 
the  form  of  a  (small)  set  of  parameters,  while  the  FT  yields  as  many  points  as  the 
original  data  set.  Taxget  identification  techniques  generally  require  a  small  number 
of  data  elements  in  order  to  be  computationally  feasible  [18].  ARMA  techniques 
supply  a  small  number  of  parameters  (the  2na  p,  and  d*  coefficients).  FT  methods 
are  nonpaxametric,  and  the  number  of  output  data  points  produced  is  the  same  as 
the  number  of  frequency  measurements  ( N ). 


3.  Data  Length  and  Decimation  Effects 


This  chapter  considers  the  effects  of  the  number  of  data  points  (frequency  sam¬ 
ples)  N  and  of  the  frequency  separation  6/  of  these  samples  on  the  estimation 
results. 

3.1  Finite  Data  Effects 

First  we  consider  the  effect  of  the  number  of  data  points,  i.e.  in  the  following 
discussion  we  assume  6j  to  be  a  constant. 

Both  the  ARMA  methods  and  the  FT  methods  are  theoretically  formulated  for 
the  infinite  data  case,  which  means  that  they  have  to  be  adapted  for  the  case  of  a 
finite  number  of  data  points.  One  way  of  looking  at  these  adaptations  is  as  follows: 
given  a  (finite)  sequence  {y(/*)}*lo\  find  an  infinite  sequence  {y(/t)}£L0  such  that 
y(fk)  —  y{fk )  ■  VfceO ...  AT  —  1,  and  such  that  the  finite  transform  (FT  or  ARMA) 
of  y  equals  the  infinite  transform  of  y. 

For  the  FT  methods  this  is  very  simple:  one  chooses  y(fk)  =  0  for  all  k  <  0 
or  k  >  N  [7];  that  is  to  say,  a  finite  FT  implicitly  assumes  that  all  measurements 
outside  the  given  data  range  are  identically  zero.  It  is  physically  obvious  that  this 
is  not  a  very  good  assumption. 

The  ARMA  methods,  on  the  other  hand,  assume  that  the  data  outside  the  known 
range  is  a  continuation  of  the  known  data  [17].  In  particular,  the  extension  y(fk) 
is  given  by  equation  (2.18)  for  k  <  0  or  k  >  N.  This  assumption  gives  a  smooth 
extension  of  the  frequency  data,  and  is  much  more  in  keeping  with  the  physical 
properties  of  the  problem  than  is  the  FT  assumption.  However,  this  assumption 
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is  valid  only  if  the  bandwidth  of  the  radar  measurements  is  small  enough  for  the 
measurements  to  be  stationary. 

A  second  way  of  looking  at  the  finite  methods  is  to  consider  the  given  sequence 
{y(fk)}^=o  as  a  subset  of  a  (unknown)  sequence  {y{fk}T=  o  an<^  describe  how  the  infi¬ 
nite  sequence  y  must  be  manipulated  to  get  the  sequence  y  described  above.  For  the 
finite  FT  methods  this  manipulation  consists  of  multiplication  with  a  rectangular 


window 


{“>*}”=  o  • 


Wj t  =  1  ,  k  —  0 . . .  N  —  1 
Wk  =  0  ,  otherwise 


This  means  that  estimated  impulse  response  will  be  the  actual  impulse  response 
convolved  with  the  transform  W(z)  of  the  window  w^. 

N-i  1  _  -N 

W(z)  =£*-*  =  — — : r  (3.1) 

k—0  1  Z 

Evaluated  on  the  unit  circle,  W(z)  is  a  sine  function,  so  the  convolution  will  have 
the  effect  of  smearing  the  estimate:  this  accounts  for  the  resolution  limit  of  the 
finite  FT  methods.  For  proofs  and  details  about  the  various  finite  FT  properties 
referred  to  above,  see  [7,8]. 

The  ARMA  method  does  not  assume  that  the  data  points  outside  the  known 
range  are  zero;  instead,  a  smooth  continuation  of  the  known  data  is  assumed.  As 
a  result,  there  is  no  windowing  effect,  and  hence  no  corresponding  resolution  limit 
for  the  ARMA  method.  Of  course,  with  both  ARMA  and  finite  FT  methods,  any 
information  contained  in  the  unknown  data  will  not  be  taken  into  account,  so  the 
ARMA  estimates  will  improve  in  quality  as  the  amount  of  available  data  increases. 

For  the  FT  methods,  the  statistical  behaviour  of  the  estimates  has  been  well- 
studied,  and  various  methods  for  trading  off  bias  for  variance  have  been  developed 


[14].  For  the  ARM  A  method,  the  estimated  parameters  in  the  model  become  ran¬ 
dom  variables.  The  statistical  behaviour  of  parameter  estimates  have  also  been 
studied  (see  [19,20]).  One  important  result  is  that  the  variance  of  the  estimates  de¬ 
creases  as  the  number  of  data  points  increases;  this  is  another  argument  supporting 
the  use  of  more  data  points  in  the  estimates. 


3.2  Decimation  Effects 


In  this  section  we  discuss  decimation  effects;  here  the  value  of  N8j  is  assumed 
constant  throughout  the  following  discussion. 

If  we  decimate  our  data,  i.e.  increase  6j  (and  therefore  decrease  TV),  we  decrease 
the  maximum  unambiguous  range  R,  since  R  =  c/(26f)  (2.14).  For  the  finite  FT 
method  this  does  not  change  the  resolution  in  meters,  which  is  determined  only 
by  the  number  of  data  points.  For  the  ARMA  methods  there  is  an  advantage 
to  choosing  the  maximum  unambiguous  range  as  small  as  possible:  if  R  is  much 
larger  than  the  size  of  the  target,  the  impulse  response  will  look  like  a  number  of 
closely  spaced  peaks,  surrounded  by  stretches  of  zero  signal  (“empty  space”).  Since 
zero  intervals  are  not  well  represented  by  a  sum  of  complex  exponentials,  ARMA 
methods  will  work  best  when  the  target  fills  the  unambiguous  range.  Therefore, 
decimation  techniques  (which  decrease  the  unambiguous  range)  are  of  interest. 

The  simplest  method  of  decimating  is  to  use  only  one  out  of  every  n  data  points, 
where  n  is  the  decimation  factor.  Of  course,  this  has  the  disadvantage  that  only 
a  fraction  of  the  available  data  is  used.  An  alternative  method  is  to  compute  the 
autocorrelation  sequence  (equation  (2.20))  using  all  available  data,  and  then  deci¬ 
mate  the  autocorrelation  sequence,  instead  of  the  data  sequence.  This  method  has 
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the  advantage  of  using  all  data,  thus  producing  more  accurate  parameter  estimates. 
This  method  is  equivalent  to  averaging  n  adjacent  data  points  to  produce  one  deci¬ 
mated  point.  In  other  words,  this  is  convolution  by  a  rectangular  window,  followed 
by  decimation.  From  the  signal  processing  literature  [8]  it  is  well  known  that  the 
rectangular  window,  while  computationally  simple,  produces  strong  sidelobes  in  the 
resulting  estimate.  Other,  more  complicated  windows  produce  far  weaker  sidelobes, 
and  thus  give  better  qua,rty  estimates  from  the  same  data.  Therefore  we  have  de¬ 
cided  to  use  a  Hamming  window  in  all  our  estimates.  Using  windows  has  the  added 
advantage  that  any  spurious  effects  in  the  data  will  tend  to  be  averaged  out. 


4.  Application  of  the  ARMA  Techniques  to  the  ESL  Compact  Range 

Data 

In  this  chapter  we  describe  the  application  of  the  techniques  developed  above 
to  the  OSU  ESL  Compact  Range  radar  data.  As  described  in  Section  2.1,  we  have 
radar  measurements  in  the  5-92  MHz  range  for  five  commercial  aircraft,  namely  the 
Boeing  707,  727  and  747,  the  DC10  and  the  Concorde.  This  data  has  been  obtained 
by  measuring  the  response  of  scaled  models  of  these  aircraft;  the  measurement 
frequencies  span  the  1-12  GHz  range. 

4.1  Explanation  of  the  Figures 

This  section  describes  the  figures  referred  to  in  the  next  sections;  it  discusses 
the  meaning  of  the  various  quantities  shown  and  the  meaning  of  the  scales  on  the 
axes. 

The  figures  that  display  actual  data  as  it  is  stored  in  the  database  are  in  units 
of  Volts/meter  (in  other  words  signal  strength,  not  power)  on  the  vertical  axis,  with 
the  index  into  the  database  file  on  the  horizontal  axis.  The  frequency  corresponding 
to  index  k  is  fk+i- 

The  plots  that  give  FFT  results  direct  from  the  data  are  in  units  of  the  signal 
power  (in  dB  meter2)  on  the  vertical  axis,  with  the  FFT  index  on  the  horizontal  axis 
The  FFT  index  k  corresponds  to  range  jf^R.  All  FFT  are  obtained  by  padding 
the  data  with  zeroes  to  obtain  256  points. 

The  other  plots  all  axe  in  units  of  power  (in  dB  meter2)  on  the  vertical  axis,  and 
range  (in  meters)  on  the  horizontal  axis.  The  range  is  corrected  for  the  two-way 


propagation  effect,  and  has  been  plotted  such  that  the  physical  center  of  the  target 
is  at  half  range;  that  is,  if  the  horizontal  scale  runs  from  zero  to  72  meters,  then 
the  center  of  the  target  is  at  36  meters.  The  center  of  the  target  is  the  point  on 
the  model  where  it  was  attached  to  the  pedestal  in  the  Compact  Range  facility; 
this  mounting  point  has  been  chosen  to  be  approximately  the  target’s  geometrical 
center.  All  ranges  have  been  corrected  for  the  target  scale  factor,  that  is  the  ranges 
given  are  in  meters  for  the  actual  aircraft,  not  for  the  model. 

For  the  ARMA  results  we  show  two  plots  per  figure,  namely  the  estimated  re¬ 
sponse  and  the  estimated  scattering  centers  (labelled  “response”  and  “scatterers”, 
respectively).  Both  are  to  the  same  horizontal  scale  to  make  comparisons  easy.  The 
“response”  plot  is  the  impulse  response  magnitude,  and  compares  to  the  FFT.  The 
“scatterers”  plot  is  a  graphical  presentation  of  the  pole  and  residue  coefficients:  each 
horizontal  line  represents  a  scattering  center;  the  vertical  tick  mark  gives  the  esti¬ 
mated  range  of  the  scattering  center,  and  the  width  of  the  horizontal  line  represents 
the  3dB  spatial  dispersion  of  the  scattering  center  (analogous  to  3  dB  bandwidth). 
The  height  of  each  horizontal  line  gives  the  energy  associated  with  the  scattering 
center. 

4.2  Preliminary  Modeling  Results 

In  this  section  we  study  a  few  preliminary  modeling  and  estimation  issues.  In 
particular,  we  study  the  effect  of  the  frequency  range  used  on  the  quality  of  the 
resulting  estimates,  the  effect  of  model  order,  and  we  verify  that  the  various  data 
preprocessing  strategies  do  not  cause  undesirable  side  effects  on  the  results.  In  ad¬ 
dition,  we  compare  the  ARMA  modeling  results  with  FT  estimates  for  verification. 
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All  results  shown  in  this  section  were  obtained  with  vertical  polarization  trans¬ 
mitted  and  received  (VV  data),  and  using  targets  at  an  aspect  angle  of  10°. 

4.2.1  Data  Decimation  Effects 

To  see  the  effect  of  using  all  available  data,  we  need  to  calculate  the  maximum 
unambiguous  range  R.  To  do  this  we  need  6/,  which  can  be  found  by  dividing 
the  Sj  of  the  database  by  the  scale  factor  of  the  aircraft  under  consideration;  for 
the  case  of  the  Boeing  707  data  this  gives  Sj  =  50  •  106/150  =  .33  •  106  Hz.  From 
equation  (2.14),  with  6/  =  0.33  MHz  we  have  that  using  all  available  data  gives  an 
unambiguous  range  of  about  450  m,  see  Figs  A.l  and  A. 2.  Both  the  FT  and  the 
ARMA  result  clearly  show  the  presence  of  the  target  (Boeing  707)  in  the  middle  of 
the  range,  and  no  response  to  either  side  of  it.  To  improve  resolution,  we  need  to 
decimate  our  data  until  we  have  decreased  the  unambiguous  range  to  approximately 
equal  the  length  of  the  aircraft.  Since  the  aircraft  is  about  70  m  long,  decimation 
of  450/70  as  6.49  is  needed;  thus,  we  need  to  use  about  211/6.49  =  32  data  points 
for  full  bandwidth  cases  (using  all  frequency  data),  or  about  16  for  half  bandwidth 
cases  (using  either  the  upper  or  the  lower  half  of  the  measurements). 

4.2.2  Hamming  Window  Data  Averaging 

The  next  set  of  figures  shows  the  effects  of  averaging  the  data  with  a  three-point 
Hamming  window.  Hamming  window  averaging  is  equivalent  to  multiplying  the 
target  response  (in  the  range  domain)  with  the  inverse  FT  of  a  Hamming  window, 
it  is  readily  verified  that  a  three  point  Hamming  window  has  a  main  lobe  width 
which  is  wider  than  any  target  of  interest.  Thus,  this  averaging  has  the  effect  of 
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attenuating  any  spurious  response  outside  the  target  range  [2]. 

Fig  A. 3  is  the  raw  data  (amplitude  and  phase)  for  one  of  the  targets  (code  name 
target  32)  The  horizontal  axis  shows  the  value  of  the  index  into  the  data  array, 
the  lowest  index  corresponds  to  1  GHz,  the  highest  index  to  12  GHz  (these  are  the 
scaled  data  base  frequencies,  not  the  frequencies  for  a  full-size  target).  Fig.  A. 4 
shows  the  windowed  data  (three  point  Hamming  window)  in  the  same  format.  As 
can  be  seen  by  comparing  the  two  Figures,  the  Hamming  window  does  not  greatly 
alter  the  data.  Figure  A. 5  shows  the  decimated  data  (decimated  by  a  factor  of  6). 
The  frequency  range  is  the  same;  i.e.,  index  1  corresponds  to  1  GHz,  and  index  36  to 
12  GHz.  The  values  in  this  plot  are  obtained  through  the  same  Hamming  window 
as  used  in  Fig.  A.4;  the  effect  of  the  decimation  can  be  readily  seen  by  comparing 
the  two  Figures:  Figure  A. 5  has  the  same  overall  shape  as  figure  A.4,  but  does  not 
have  the  more  rapid  variations  seen  in  the  latter  figure. 

Figure  A. 6  shows  the  result  of  an  FT  on  the  raw  data  of  Fig.  A. 3.  This  clearly 
shows  the  target  (the  bumps  at  the  edges  of  the  plot),  and  the  expected  empty 
space  around  it  (the  low  region  in  the  middle).  Figure  A. 7  shows  the  result  of  an 
FT  on  the  windowed  data  of  Figure  A.4.  By  comparing  with  the  FT  result  for 
the  raw  data  (Fig.  A. 6),  it  is  clear  that  the  window  achieves  its  intended  function 
of  reducing  the  spurious  responses  in  “empty  space”  region  without  altering  the 
response  in  the  “aircraft”  region:  this  will  reduce  aliasing  effects  when  the  data 
is  decimated.  Finally,  Figure  A. 8  shows  the  FT  of  Figure  A. 5,  the  windowed  and 
decimated  data.  We  can  clearly  see  that  the  “aircraft”  region  has  now  expanded  to 
occupy  the  entire  plot  range,  as  intended. 


4.2.3  Stationarity  of  the  Data 


Going  back  to  the  previous  example  (the  Boeing  707),  Figure  A. 9  shows  the 
ARMA  results  when  we  use  frequency  data  in  the  6.67-40  MHz  range,  about  half  of 
the  available  bandwidth;  Figure  A. 10  is  the  corresponding  FT.  Figure  A. 11  shows 
the  ARMA,  and  Figure  A. 12  shows  the  FT  results  for  data  from  the  40-  80  MHz 
range,  the  upper  half  of  the  frequency  band.  Comparison  between  the  ARMA 
response  curves  and  the  FT  curves  shows  that  they  agree  only  marginally  well.  One 
reason  for  this  is  that  the  FT  curve  represents  the  ideal  response  convolved  with  a 
window  function;  the  ARMA  response  does  not  contain  this  windowing  operation. 
Second,  it  is  the  area  under  a  peak  in  the  ARMA  curve  that  represents  its  energy; 
thus,  sharp  peaks  with  little  energy  seem  to  disagree  with  FT  results,  when  in 
fact  the  agreement  is  far  better.  We  have  found  it  more  helpful  to  compare  the 
scatterers  plots  with  the  FT  plots;  here  one  sees  direct  correspondence  between 
scattering  centers  of  significant  energy  and  peaks  in  the  FT  curves. 

To  consider  “data  stationarity”,  compare  the  lower  frequency  plots  with  the 
higher  frequency  plots.  It  can  be  seen  that  the  ARMA  estimates  for  the  lower 
frequency  range  give  much  more  spatially  distributed  scatterers  (that  is  scatterers 
that  are  not  point-like,  but  occupy  a  non-zero  range).  These  results  imply  that: 

•  the  data  is  not  stationary  with  respect  to  frequency.  This  is  to  be  expected 
from  electromagnetic  theory,  but  we  have  now  observed  that  the  effect  is 
significant  in  the  frequency  ranges  and  bandwidths  used  in  this  study. 

•  Since  the  ARMA  model  is  well  suited  to  estimating  point-like  scatterers  (i.e. 
a  more  optical  response),  we  should  use  higher  frequency  data  whenever  pos- 


sible,  because  more  point-like  scatter-behaviour  is  seen  in  this  range.  It  is 
expected  that  data  taken  at  higher  frequencies  will  result  in  better  perfor¬ 
mance  of  the  ARMA  techniques. 

Based  on  these  observations,  all  remaining  experiments  use  data  in  the  40-80  MHz 
range. 

4.2.4  Removal  of  Spatially  Distributed  Scattering  Centers 

To  further  reduce  the  effects  of  non-stationarity,  and  improve  performance,  we 
use  a  technique  that  has  the  added  advantage  of  enabling  the  ARMA  estimator  to 
automatically  choose  the  number  of  scattering  centers  to  include  in  the  model;  for 
estimation  on  unknown  target  data  this  is  an  important  and  desirable  property  of 
the  algorithm. 

The  technique  uses  the  fact  that  the  estimation  procedure  entails  two  main  steps: 
first  the  locations  and  spatial  distribution  of  the  scattering  centers  are  estimated, 
then  their  energies  are  computed.  The  modified  algorithm  performs  the  first  step 
with  a  fixed  number  na  of  scattering  centers.  Then  all  centers  that  are  too  spatially 
distributed  are  eliminated.  Finally,  the  amplitudes  of  the  remaining  centers  are 
estimated  as  before. 

The  spatial  distribution  of  a  scattering  center  is  related  to  the  magnitude  of  the 
pole  estimated  for  that  center  (see  (2.17));  the  contribution  of  the  center  to  the 
frequency  response  at  frequency  /,  is  (C.f.  (2.18)) 

dtPi  0  <  k  <  N  —  l 

If  |p,|  is  too  far  from  unity  (i.e.  if  the  scattering  center  is  too  spatially  distributed), 
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then  its  contribution  to  the  frequency  model  is  greatly  different  at  frequency  /0  then 
at  frequency  /jv-i-  The  ratio  of  these  contributions  is: 


(4.1) 


d  =  V  M"-'  .  IPil>l 

d  =  -£r  ip,rw+l  ,  m  <  i 

In  the  modified  algorithm,  if  |p,  |  is  such  that  D  is  greater  than  some  threshold,  that 
pole  is  discarded. 

Empirical  tests  showed  that  D  =  100  provided  good  removal  of  highly  dispersive 
scatterers  without  removing  sharp  scatterers.  As  an  example,  Figures  A. 13  and  A. 14 
show  the  results  for  D  =  1000  and  D  =  100,  respectively.  It  can  be  seen  that  the 
results  are  very  similar  for  the  two  cases.  Also,  comparing  Figure  A.  14  with  Figure 
A.  11  (where  no  pole  removal  is  used)  shows  that  sharp  scatterers  are  more  clearly 
seen  if  D  =  100  is  used.  Finally,  it  was  found  that  estimation  results  were  less 
sensitive  to  model  order  selection  when  this  technique  is  used,  as  discussed  below. 


4.2.5  Effects  of  Model  Order 

The  ARMA  estimator  is  based  on  the  assumption  that  the  target  can  be  modeled 
as  a  set  of  na  scattering  centers,  where  the  value  of  na  needs  to  be  given.  This  section 
discusses  the  effect  different  choices  of  na  has  on  the  ARMA  estimates. 

Figures  A. 15  through  A. 21  show  the  effect  of  increasing  the  order  of  the  ARMA 
estimate  from  1  to  9;  compare  with  the  FT  response  in  Figure  A. 12.  The  pole 
elimination  procedure  was  used  in  these  estimates  with  D  =  100,  so  the  number  of 
scatterers  in  thew  estimate  is  not  equal  to  the  original  model  order  in  many  cases. 
It  should  be  noted  that  for  orders  1  and  3  the  estimate  did  not  contain  any  useful 
scatterers  (all  scaiterers  were  too  distributed)  therefore,  no  plots  are  included  for 


these  two  cases.  From  the  Figures  we  can  clearly  see  that,  as  early  as  at  fourth  order 
the  estimater  picks  out  the  two  major  peaks  in  the  FT,  which  appear  to  correspond 
to  the  cockpit  windows  and  the  tail  assembly.  At  sixth  and  higher  orders  we  see 
additional  peaks  which  seem  to  correspond  to  the  leading  edge  of  the  engines  and 
the  trailing  edges  of  the  wing  and  the  engines.  At  some  orders  (order  4,  for  example) 
we  get  a  peak  that  seems  to  be  a  post-resonance. 

The  correspondence  between  the  responses  and  the  physical  atributes  of  the 
target  can  be  seen  by  overlaying  Figure  A. 113  (a  scaled  picture  of  the  Boeing  707), 
at  an  aspect  angle  of  10°,  on  the  plots. 

In  evaluating  these  results,  we  should  keep  in  mind  that  the  Boeing  707  is  a 
geometrically  complicated  target:  it  has  four  engines,  a  complex  tail  assembly, 
complex  wing  shape,  and  several  irregularities  on  the  body  (cockpit  windows,  an¬ 
tennas).  Such  a  complex  target  needs  a  high  model  order  (for  good  modelling,  we 
need  the  model  order  to  be  at  least  as  high  as  the  number  of  features,  or  about  15 
for  this  plane),  but,  due  to  the  small  number  of  data  points,  we  cannot  increase 
the  model  order  that  far:  in  other  words,  this  is  a  hard  case  for  the  ARMA  estima¬ 
tor.  Despite  this,  the  ARMA  estimator  seems  to  give  reasonable  results,  and  the 
automatic  elimination  of  spatially  distributed  scatterers  seems  to  work  well. 

4.3  Aircraft  Data  Results 

In  this  section  we  present  two  sets  of  results.  The  first  set  shows  two  aircraft, 

the  Concorde  and  the  Boeing  707,  rotated  through  aspect  angles  0° - 180°,  both 

at  horizontal  transmitted  and  horizontal  received  polarization  (HH).  The  second  set 
shows  all  five  aircraft  at  an  aspect  angle  of  10°,  at  both  HH  and  VV  polarization. 


4.3.1  Comparison  of  Estimates  at  Different  Aspect  Angles 

We  expect  from  geometrical  considerations  that  as  an  aircraft  rotates,  its  re¬ 
sponse  will  vary  more  or  less  smoothly  as  a  function  of  angle,  so  the  aircraft  will 
look  different  at  different  angles,  but  similar  at  similar  angles.  The  DFT  results  for 
the  Concorde  (Figures  A. 79  through  A. 97)  do  indeed  exhibit  this  behaviour,  the 
DFT’s  for  the  Boeing  707  (Figures  A.41  through  A. 59)  do  this  to  a  lesser  degree. 
The  ARMA  results  for  the  Concorde  (Figures  A. 60. .A. 78)  show  this  behaviour  for 
angles  close  to  0°  and  180°,  but  not  at  intermediate  angles.  For  example,  the  peak 
at  range  48  m  in  Figure  A. 60  can  be  seen  moving  smoothly  toward  40  m  as  the 
angle  varies  from  0°  to  30°.  Except  for  a  few  angles,  the  ARMA  results  for  the 
Boeing  707  do  not  show  this  behaviour  at  all.  We  suspect  this  is  caused  by  the  fact 
that  the  ARMA  model  attempts  to  extract  the  most  characteristic  features  of  the 
response,  rather  than  model  it  in  its  entirety,  and  the  relative  importance  of  various 
scatterers  does  not  have  to  vary  smoothly  with  the  aspect  angle. 

Despite  its  non-continuous  behaviour,  the  ARMA  model  seems  to  follow  the 
most  important  peaks  of  the  DFT  rather  faithfully,  with  exceptions  at  110°  and 
140°  for  the  Concorde  and  at  0°,  20°,  50°,  80°,  90°  and  140°  for  the  Boeing  707. 
Again  it  should  be  observed  that  the  Boeing  707  is  a  complex  plane,  which  is 
somewhat  difficult  to  model. 

4.3.2  Comparison  of  Estimates  from  Different  Aircraft 

The  second  set  of  results  shows  all  five  planes  (Boeing  707,  727,  747;  DC10 
and  Concorde)  at  10°  aspect  (10°  is  chosen  because  0°  is  a  degenerate  case  where 
the  complexity  of  the  target  is  reduced  through  symmetry)  in  both  Horizontal  and 
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Vertical  polarization. 

As  can  be  seen  by  comparing  Figures  A. 103  through  A. 107  and  A. 108  through 
A. Ill  (the  ARMA  results)  with  Figures  A.98  through  A. 102  (the  DFT  results),  we 
again  have  a  fairly  good  correspondence  between  the  DFT  and  the  ARMA  results. 
It  is  also  easy  to  see  that  the  ARMA  parameters  for  different  aircraft  are  different, 
so  a  classifier  based  on  ARMA  results  should  be  able  to  distinguish  between  these 
planes.  Quantative  analysis  of  such  an  ARMA-based  classification  procedure  is  a 
topic  of  future  research. 
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5.  Conclusions 


Prom  our  simulations  and  analysis,  we  can  conclude  the  following: 

1.  ARM  A  models  can  be  used  to  estimate  the  impulse  response  of  a  radar  tar¬ 
get,  given  radar  cross-section  measurements  at  a  number  of  different  (evenly 
spaced)  frequencies.  Our  simulations  show  that  about  20  frequencies  is  enough 
to  get  satisfactory  results. 

2.  The  resulting  ARMA  estimates  have  the  advantages  over  estimates  obtained 
through  Fourier  Transform  techniques  that 

•  They  tend  to  concentrate  on  the  strongest,  most  pointlike  scattering 
centers  in  the  response. 

•  They  are  not  resolution-limited  by  the  bandwidth  of  the  measurements. 

•  They  are  expressed  in  the  form  of  a  small  number  of  parameters,  rather 
than  as  many  samples  of  a  curve. 

Spacing  of  the  frequency  samples  should  be  chosen  such  that  the  maximum 
unambiguous  range  of  the  target  response  is  at  least  as  large  as  the  target.  If  the 
spacing  of  the  frequency  samples  is  chosen  such  that  the  maximum  unambiguous 
range  is  much  larger  than  the  target  length,  modelling  quality  can  be  improved 
by  applying  data  decimation  techniques  before  the  actual  estimation.  In  a  noisy 
environment,  such  a  strategy  could  be  used  to  improve  the  signal  to  noise  ratio. 

The  total  bandwi  th  of  the  measurements  should  be  restricted  so  as  to  avoid 
violating  the  implicit  stationarity  assumption  of  the  ARMA  model.  For  the  cases 
we  have  studied,  we  found  that  a  bandwidth  of  40-80  MHz  works  well. 
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While  no  formal  qualitative  studies  have  been  made,  our  simulations  indicate 
that  the  ARMA  modelling  technique  is  capable  of  distinguishing  between  different 
aircraft,  and  between  different  aspect  angles  for  the  same  aircraft.  In  other  words, 
the  ARMA  modelling  technique  seems  to  be  a  suitable  signal  processing  method 
to  precede  target  classification.  Future  research  should  focus  on  a  more  formal, 
quantative  evaluation  of  classification  performance  using  ARMA  parameters. 
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A.  Figures 


The  following  pages  contain  the  Figures  for  this  report 
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Figure  A.l:  ARMA  Response  of  Boeing  707  at  Vertical  Polarization,  10°  aspect 
angle,  40-80  MHz,  120  frequency  samples  used 
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Figure  A. 6:  FT  of  Target  number  32,  Vertical  Polarization,  10°  aspect  angle 
l-12GHz  model  frequencies,  all  frequency  samples  used. 
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Figure  A. 7:  FT  of  Target  number  32,  Vertical  Polarization,  10°  aspect  angle, 
l-12GHz  model  frequencies,  with  3-point  Hamming  window,  all  frequency  samples 
used. 
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Figure  A. 9:  ARMA  Response  of  the  Boeing  707  at  Vertical  Polarization,  10°  aspect 
angle,  16  frequency  samples  used,  6.67-40  MHz. 


41 


Figure  A.  13:  ARM  A  Response  of  the  Boeing  707  at  Vertical  Polarization,  10°  aspect 
angle,  20  frequency  samples  used,  40-80  MHz,  rejection  factor=1000 


Figure  A.  14:  ARM  A  Response  of  the  Boeing  707  at  Vertical  Polarization,  10°  aspect 
angle,  20  frequency  samples  used,  40-80  MHz,  rejection  factor=100 


Figure  A.  15:  ARMA  Response  of  the  Boeing  707  at  Vertical  Polarization,  10°  aspect 
angle,  20  frequency  samples  used,  40-80  MHz,  2nd  order  model 


Figure  A. 23:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  10° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 24:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  20° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A.31:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  90 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 


Figure  A. 33:  ARM  A  Response  of  the  Boeing  707  at  Horizontal  Polarization,  110 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 


Figure  A. 34:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  120 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 


Figure  A. 37:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  150 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 39:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  170° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 40:  ARMA  Response  of  the  Boeing  707  at  Horizontal  Polarization,  180° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 41:  FT  of  the  Boeing  707  at  Horizontal  Polarization,  0°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A. 43:  FT  of  the  Boeing  707  at  Horizontal  Polarization,  20 
frequency  samples  used,  40-80  MHz 


Figure  A.44:  FT  of  the  Boeing  707  at  Horizontal  Polarization,  30°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A. 45:  FT  of  the  Boeing  707  at  Horizontal  Polarization,  40°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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re  a.di:  M  ol  the  Boeing  707  at  Horizontal  Polarization,  100°  aspect  angle, 
equency  samples  used,  40-80  MHz 
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Figure  A. 60:  ARMA  Response  of  the  Concorde  at  Horizontal  Polarization,  0°  aspect 
angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 63:  ARM  A  Response  of  the  Concorde  at  Horizontal  Polarization,  30° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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|  Figure  A. 64:  ARMA  Response  of  the  Concorde  at  Horizontal  Polarization,  40° 

|  aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 65:  ARMA  Response  of  the  Concorde  at  Horizontal  Polarization,  50 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 73:  ARMA  Response  of 
aspect  angle,  20  frequency  samples 
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Figure  A. 76:  ARMA  Response  of  the  Concorde  at  Horizontal  Polarization,  160° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A. 79:  FT  of  the  Concorde  at  Horizontal  Polarization,  0°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A. 83:  FT  of  the  Concorde  at  Horizontal  Polarization,  40°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A. 86:  FT  of  the  Concorde  at  Horizontal  Polarization,  70' 
frequency  samples  used,  40-80  MHz 
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Figure  A.89:  FT  of  the  Concorde  at  Horizontal  Polarization,  100°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A.91:  FT  of  the  Concorde  at  Horizontal  Polarization,  120°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A. 95:  FT  of  the  Concorde  at  Horizontal  Polarization,  160°  aspect  angle,  20 
frequency  samples  used,  40-80  MHz 
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Figure  A. 100:  FT  of  the  Boeing  747  at  Horizontal  Polarization,  10°  aspect  angle, 
20  frequency  samples  used,  40-80  MHz 
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Figure  A.  105:  ARM  A  Response  of  the  Boeing  747  at  Horizontal  Polarization,  10° 
aspect  angle,  20  frequency  samples  used,  40-80  MHz,  9th  order  model 
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Figure  A.  113:  The  Boeing  707 
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